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Abstract: 


This paper introduces a new approach aimed at restoring images corrupted by random valued 
impulse noise. The adopted methodology leverages fuzzy logic and encompasses three primary 
stages: estimation of noise density, detection of fuzzy noise, and reduction of fuzzy noise. Within 
the fuzzy noise detection phase, a fuzzy set labeled as "Noise-Free" is formulated through the 
utilization of the rank-ordered mean of absolute differences and the estimated noise density. This 
set serves to discern whether a given pixel should be classified as noisy or noise-free. Utilizing the 
fuzzy logic in the proposed method collaborates to determine the ultimate fuzzy weight assigned to 
each pixel, thereby facilitating the restoration of corrupted image pixels. Empirical results based on 
peak signal-to-noise ratio, mean square error, and visual assessment demonstrate the effectiveness 
of the proposed technique in suppressing noise, preserving fine details, and surpassing the 
performance of several established filtering methods. 


1. Introduction 


In numerous digital image processing applications, the captured images often suffer from noise 
corruption, stemming from issues during image transmission or acquisition. Noise degradation 
diminishes image quality and introduces adverse effects into subsequent image processing stages, 
including segmentation, parameter estimation, and enhancement. Consequently, the elimination of 
noise from images emerges as a fundamental and critical task within the domain of image 
processing. Impulse noise, characterized by brief 'on/off’ noise pulses of relatively short duration, 
holds a significant and pervasive presence among digital image noise sources. Its impact is felt 
during image acquisition due to factors like noisy sensors (arising from switching or sensor 
temperature variations), as well as during transmission due to channel imperfections (interference, 
atmospheric disturbances), hardware-related defects, or synchronization errors during analog-to- 
digital conversion within image processing workflows. This study specifically focuses on the 
Random Valued Impulse Noise (RVIN) model of impulse noise [1-5]. 


Among the diverse techniques for mitigating impulse noise, the median filter [6,7] has gained 
popularity due to its effective noise suppression capabilities coupled with high computational 
efficiency. However, the median filter often suffers from drawbacks such as destroying image 
details and introducing blurring, as each pixel is replaced with the median value of its surrounding 
neighborhood. To address this issue, several forms of adaptive and center-weighted median filter 
[8-15] was introduced, which assigns greater emphasis to the central pixel within the observed 
window. Although this modification preserves more image details than the median filter, it 
uniformly treats all pixels across the image, without distinguishing between noisy and noise-free 
pixels. This paper introduces a new approach for removing RVIN from images. The proposed 
method employs fuzzy techniques alongside estimated noise density. The estimation of noise 
density serves dual purposes: first, to determine the appropriate window size for the fuzzy noise 
detection step, and second, to shape the membership function used in the fuzzy noise detection 


process. The fuzzy logic, which is adopted by Zadeh in 1965, is a mathematical tool for dealing 
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with uncertainty. It presents a crucial model of computing with words and furnishes an approach 
for handling imprecision and information granularity within the realm of soft computing [16-20]. 

2. RVIN Model 

Impulse noise inherently maintains its independence and lacks correlation with individual image 
pixels, resulting in a subset of pixels bearing noise while the remainder remains untarnished. When 
dealing with random valued impulse noise, colloquially termed uniform impulse noise, the affected 
pixels can adopt any value within the gray level spectrum, denoted as [Lmin, Lmax] in an 8-bit 
image representation. Consequently, these noisy pixels span the full gray level range from 0 to 255. 
This scenario entails a stochastic dispersal of noise throughout the entire image, leading to an equal 


probability of any given gray level value manifesting as noise. The Random Valued Impulse Noise 
(RVIN) model can be briefly described as follows [21, 22]: 


n(x, y), with probability P (1) 
1G y) = arr 
O(x, y), with probability 1 — P 


Where, (x,y) is the noisy image pixel, n(x,y) is the noisy impulsive pixel at position (x,y), and 
(x,y) is the uncorrupted (original) image pixel. 
3. Proposed Method 


The proposed method consists of four main concepts: estimation of noise density, noise detection, 
noise removing, and iterative processing which are described below in details. 


3.1. Estimation of Noise Density 


For estimating RVIN density in color image, each channel (color component) is processed 
separately by dividing it into 16 non-overlapping blocks (Bi) with i € {1,2,...,16} and then the 
noise density for each block in each channel is calculated as described below for the first block 
(B,) only but it is performed in an analogous way for the others: 


1- Calculating the standard deviation (std,) for the entire block B, based on the following 
equation: 


std, = laa D ves ha y)- m,)- we (2) 


where, ™m, : is the mean ane of first block By; 
I(x, y) : is the pixel data inside B, . 


2- Scanning the whole image block B, using a 3x3 sliding window from pixel to pixel and using 
the following three steps in each scan: 


i) DN =|i(x+s,y +t) —I(x,y)| withs,t € Syy (3) 


Where, Dy represents the set of absolute differences between the central pixel and its four nearest 
neighbors within the processed window and 
Sy4 denotes the set of coordinates for the four neighbors of I(x, y) given by: 


Sn = {(-1,0), (1,0), (0, =i), (0,1)} see (4) 
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Then the elements in Dy can be arranged in ascending order by value: 


1 2 3 4 
24 20.20, 


ii) Defining v,yy which represents the mean value of the second and third ordered elements in 
D}y¥ by the following equation: 
_ dhyt+dRy 


Vgy = OY 30) 


iii) Using the standard deviation (std) of the processed image block and the value of vy,y to 
identify the current pixel I(x, y) as follows: 


1, Yy% > std, 


Mpy (X,Y) = .. (6) 
0, Vyy S std, 


Where, Mpy (x,y) denotes the decision rule for estimating RVIN density. 


3- The noise density NDpy for the first block is given as: 


d«y)eB, Mav(~%y) 


NDB. = CF x x x 100% oe, 


By 


Where Mg, X Ng, denotes the size of the block B,; , and CF denotes the control factor which is 
estimated by experiments CF = 1.35 (see table 1). 


4- The calculation of the noise density for the other image blocks is done in the same way as 
for (ND fx, ) . Then, define the minimum value of (NDF, ) as: 


NDmin = min (NDB ...(8) 
am cum , with i € {1,2,...,16} 
5- If ND,»in is less than one, then the image is classified as clean. Otherwise, the image is 
classified as corrupted with RVIN. In this case, the overall noise density of the image equals to 
the mean value of the noise densities across all blocks and expressed is as follows: 


0,if NDmin <1 
NDpy = see (9) 


meany<ici6(ND ) otherwise 


Where, NDpy is the RVIN density for the entire channel. 
3.2. Nosie Detection 


As stated in section 2, in the context of RVIN, a noisy pixel holds the potential to exhibit any value 
within the gray level range [0-255] for an 8-bit image representation. Additionally, it might slightly 
deviate in intensity from the original value. Consequently, mitigating such noise presents a greater 
challenge compared to dealing with salt and pepper noise, necessitating a robust detection strategy. 
The proposed method employs the Rank-Ordered Mean of the Absolute Differences (ROMAD) in 
conjunction with estimated noise density, integrating the principles of fuzzy logic. The process of 
calculating the rank-ordered mean difference involves the central pixel's comparison with its 
neighboring pixels within a sliding window of dimensions (2K+1)x(2K+1), encompassing each 
pixel in the image. The value of K is dynamically determined based on the prevailing noise density: 
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K is set to 1 when (VN RV<35%) and otherwise, it takes on the value of 2. To compute (ROMAD) 
for a given pixel (x,y), the following three steps are executed: 


1- Compute the absolute differences, denoted as D,y between the current pixel /(x,y) and 
neighboring pixels within the observed window using the following formula : 


Duy = (x + s,y + t) —1(x, y)| with s,t € {-K,..., +K}, (s,t) # (0,0) ... ...(10) 
Where, /(x + s, y + t) denotes the pixels values within the processed window. 


2- Assume aby be the ith smallest ranked value in D,,, when the elements of D,, are arranged in 
ascending order, such that 
Cyvedy Ss Say 


Where, u is the number of the central pixel neighbors in the processed window (2K + 1) x (2K + 
1). So, u = ((2K +1)? -1) 


3- The (ROMAD) for the current pixel is defined as: 


aE 
= eed 1D) 


n 
Where, n = u — 2(2K — 1) ... (12) 


Upon completing the calculation of ROMAD for every pixel within the image, the ROMAD value 
serves as a valuable metric for discerning between pixels that are affected by noise and those that 
remain noise-free. Specifically, in scenarios where the noise ratio is low and the ROMAD value for 
a given pixel I(x,y) is small, it signifies that the pixel is noise-free and exists within a homogenous 
neighborhood. Conversely, a larger ROMAD value suggests that the pixel may be noisy or situated 
on an edge. Consequently, a fuzzy set labeled "Small ROMAD" is introduced to capture the 
linguistic notion of "small," and its membership function, denoted as sm q1 and illustrated in figure 
1, is employed to identify noise-free pixels. However, a challenge arises when dealing with high 
noise ratios, leading to consistently elevated ROMAD values. To address this issue, the 
membership function Usmaqi for the "Small ROMAD" fuzzy set can be dynamically shaped (i.e., its 
parameters determined) for each neighborhood. This adaptability is governed by two key criteria: 
firstly, the homogeneity level (H,,)of the neighborhood surrounding the tested pixel I(x,y). 
Secondly, the overall noise density of the entire image. 


R 


Usmall 
1 
“Small ROMAD” 
0 Difference 
) a b 


Fig.1: Membership function of fuzzy set “Small ROMAD” 
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The following equation describing the membership function Usmaiy : 


1,Ryy <a 
b-Ry 
Usmau(Rxy) = ==", S Rey Sb (13) 
0;Rey > b 
Where: 
a= Ayy x NFry on (14) 
b=q X NFry .. (15) 


NFpy denotes the noise factor and q represents the nth smallest value in Ra4sy+t)- 
Further, the following steps describing how H,,, , NFpy and q can be obtained: 
Firstly, H, can be computed by utilizing the ROMAD values in the observed window as: 


1- Assume 1; be the ith smallest ranked value in Ry.45y42) When the elements of Rav+s,y4t) are 
arranged in ascending order, such that: 


1nSsnsesn 
Where, u = ((2K + 1)? -1) 


2- obtain the value of H,y as: 


Axy = tenth ... (16) 
Secondly, the value of g is obtained as: 
q=T™ ... (17) 


Thirdly, the value of NFpy can be calculated by utilizing the straight line equation as illustrated in 


figure 2 and using the following expression: 
NFay = =e) _. (18) 


P3—P2 
NFry 


Py 


P, P3 WD sy 


Fig. 2: The function of NF py 


Where the parameters (P1, P2, and P3) are determined through experimental investigations. 
Empirical findings have indicated that the optimal values for these parameters are set to (2.9, 0, 61) 
respectively when the estimated noise density remains at or below 35%. Conversely, when the 
estimated noise density surpasses 35%, the parameters (P1, P2, and P3) are optimized at values of 
(1.9, 20, 100) respectively. 
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Lastly, for each tested pixel /(x,y), a fuzzy set “Noise-Free” is generated using the following 
fuzzy rule: 


Fuzzy Rule 1: Defined when a central pixel /(x, y ) is a noise free pixel: 
IF Ry, is small 
THEN the central pixel I(x, y ) is noise- free 


This rule can be applied by employing the equality operation between two fuzzy sets. As a result, 
the membership function of the fuzzy set “Noise-Free” is calculated as: 


Unoisefree (I (X, y)) = Beran Ras) sa V9) 

3.3 Noise Removing 

The noise removing step is demonstrated for the red component only (1.e., Cp) but it is implemented 
in an analogous way for the other components. Hence, assume that Cp(x,y) is a noisy pixel 
(1.€. Unoisefree(CR(x,y) <1)) and Fp(x,y) is the corresponding pixel of the filtered red 
component, and then one of the following cases will be applied for restoring the noisy pixel as 


follows: 

Casel: IF (Unoiserree(Cg(*,¥) <1) AND Unoisefree(Ca(x,y) =1)) 

Fp(x, y) = max(min(C,(x, y) + Agp(x, y), 255), 0) ... (20) 
Case2: IF (Unoisefree(Cq(*,¥) = 1) AND Unoisefree(Ca(x,y) <1)) 

Fp(x, y) = max(min (C¢(x, y) + Arg(x, y), 255) , 0) ww 21) 


Case3: IF (Pacteapree (Cox, y) =1 ) AND Unoise free (C(x, y) =1 )) 
Fa(x, y) = max(min(0.5(Cg(x, y) + Agg (x, y) + Cp(x, yy) + Ape (x, y)), 255), 0) ... (22) 
Case4: IF (Haoisapres (Cox, y) <i ) AND Hnoise free (Ca(X, y) <1 )) 


Fa(x% y) = ee 
Vs=-K Lt=—-K CR(X+5,Y+t).Unoise freeCR(X+S,y+t) 
1- F Cex, ) ak aes ve iia 
( Unoise free RC y) Woy cy Unoise freeCR(X+5,y+t) 


(apts green Os y)CR(x, y)) shes (23) 


Where, Cp(x +s, y + t) represents the pixels values in the considered window of size (2K + 1) x 
(2K + 1). The size of observed window is selected adaptively according to the number of the noise 
free pixel in that window starting with W = 1. If the observed window is fully noisy, then the size 
of window will be increased until the condition (GB > 0) is met. 

KK 


= ». > Ca(x+s,y +t) with Unoisefree(Ca(x +s,y+t)) =1 4) 
s=—-K t=-K 
3.4 Iterative Processing 


During the noise detection stage, the membership functions L,smqi, as illustrated in figure 1, 
dynamically adjust their shape based on the estimated noise density. Consequently, the proposed 
algorithms are executed in an iterative manner with two iterations. This iterative approach aims to 
enhance noise removal while preserving the details of image. In the second iteration, the altered 
image obtained from the initial iteration is utilized 
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4. Experimental Results 


Identifying a suitable parameter is a very significant task to obtain successful filtering results. 
Hence, different experiments have been implemented to select the best values for the control factor 
CF Eq. (7) which has a significant effect on the result of the noise density estimation. Table 1 
contain different noise levels with parameters varying over a range of possible values for the 
“Peppers” image. It is clear from table 1 that the best value of the parameter CF equal to 1.35, 
where the difference between the estimated noise density and the real noise density is less than 5%. 
Hence, the suggested value is CF =1.35. 


Table 1: Determination of parameter CF on “Peppers” image 


Color Parameter Real Noise Density (R VIN) 

components values 10% 20% 30% 40% 50% 60% 
0.95 7.8% | 15.0% | 22.0% | 29.2% | 35.4% | 40.1% 
Red 1.15 9.4% | 18.1% | 26.7% | 35.3% | 42.9% | 48.6% 
component 130 11.0% | 21.3% | 31.3% | 41.5% | 50.4% | 57.0% 
is) 12.7% | 24.5% | 35.9% | 47.6% | 57.8% | 65.5% 
1.75 14.3% | 27.6% | 40.6% | 53.8% | 65.3% | 73.9% 
0.95 5.9% | 12.5% | 19.4% | 25.9% | 33.0% | 38.8% 
fan 1.15 7.1% | 15.2% | 23.6% | 31.4% | 39.9% | 46.9% 
ecinpcnent 1.35 8.4% | 17.8% | 27.7% | 36.9% | 46.9% | 55.1% 
1.55 9.7% | 20.5% | 31.8% | 42.3% | 53.9% | 63.3% 
1.75 10.9% | 23.1% | 35.9% | 47.8% | 60.8% | 71.5% 
0.95 7.2% | 15.1% | 22.4% | 29.3% | 35.8% | 40.7% 
Blue 1.15 8.8% | 18.3% | 27.2% | 35.5% | 43.4% | 49.3% 
Somponeil Loo 10.3% | 21.5% | 31.9% | 41.7% | 50.9% | 57.9% 
LD 11.9% | 24.7% | 36.7% | 47.8% | 58.5% | 66.5% 
1.75 13.4% | 27.8% | 41.4% | 54.0% | 66.0% | 75.0% 


The performance evaluation of the suggested method involves a comparison with several existing 
noise reduction methods: DWM [23], HFC [24], CAFSM [25], and FMVMEF [26]. Furthermore, the 
same proposed method is applied to grayscale images without incorporating the color differences 
step, termed as proposed (gray) in this section. Table 2 presents the quantitative outcomes of 
objective quality measurements, specifically Mean Squared Error (MSE) and Peak Signal-to-Noise 
Ratio (PSNR) [27-30], for the “Peppers” image (refer to figure 4). This image is subjected to 
corruption by 10%, 30%, and 50% RVIN. It is obvious from table 1 that the proposed color image 
algorithm has a superior performance as compared with the other noise reduction methods and with 
the proposed gray scale image algorithm as well. The impact of noise density on the proposed 
algorithm's performance, as well as the performances of related methodologies in terms of PSNR, is 
visualized in figure 3. While, figure 4 illustrates the noise filtering and detail preservation of the 
proposed method as compared with other methods. It is obvious from figure 4 that the proposed 
color image algorithm obtains the best result in low and high noise densities due to the following 
main reasons: 
1. Adaptive Window Size: The utilization of an adaptive window size during both the noise 
detection and noise filtering stages is informed by the number of noise-free pixels present 
within the observed window. 
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2. Color Component Correlation: Exploiting the inter-correlation between color components 
during the noise filtering stage. So, for restoring a certain noisy pixel Cp(x, y) at position(x, y), 
the corresponding pixels in the other components i.e., (C¢(x, y) and Cg(x,y)) are used instead 
of the neighbors pixels for Cp (x, y). 


3. Robust Fuzzy Noise Detection: The integration of a robust and powerful fuzzy noise detection 
scheme further contributes to the exceptional performance of the suggested method. 


Table 2: Numerical results of the suggested method with related works using “Peppers” 


image. 
Noise density 
10% 30% 50% 

Method 
MSE MSE MSE 
PSNR | (ag) | PSNR | ag) | PSNR | Gao) 
Noisy 18.30 1.4932 13.43. | 4.5994 11.15 7.7194 
FMVMF 30.85 0.0822 | 21.53 0.7111 16.00 | 2.5806 
HFC 29.39 | 0.1036 17.75 1.6909 12.94 | 5.1493 
CAFSM 28.32 | 0.1528 26.04 | 0.2556 | 22.67 | 0.5519 
DWM 31.40 | 0.0726 | 26.65 0.2164 | 22.44 | 0.5753 
Proposes 34.81 0.0332 29.48 0.1146 25.62 0.2778 

(gray) 
Proposed 37.24 0.0191 30.82 0.0840 26.24 0.2397 

(color) 


—>— Proposed (color) 
--+-- Proposed (gray) 


Noise density (RVIN) 


Fig. 3: Comparison chart in term of PSNR of the proposed color image algorithm with 
related works using the “Peppers” image corrupted with RVIN (10% - 60%). 
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Fig. 4: Results of RVIN removing of “Peppers” image, (a) original image, (b) Noisy image 
corrupted with 30 % RVIN (PSNR: 13.43), (c) Processed image with HFC (PSNR: 17.75), (d) 
Processed image with FMVMF (PSNR: 21.53), (e) Processed image with CAFSM (PSNR: 
26.04), (f) Processed image with DWM (PSNR: 26.65) (g) Processed image with the proposed 
gray image algorithm (PSNR: 29.48), (h) Processed image using the proposed color image 
algorithm (PSNR: 30.82). 


5. Conclusion 
This study suggests a new method for RVIN detection and reduction method for color image. The 
proposed method utilizes fuzzy logic, noise density estimation, and correlation between RGB 
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components. The method capitalizes on the calculated noise density to dynamically determine an 
optimal window size for noise detection and shape the membership function employed in the noise 
detection process. Through comprehensive experimentation, the proposed method showcases 
exceptional performance surpassing numerous established impulse noise filtering techniques in 
both objective and subjective evaluations. Notably, the method excels in effectively removing noise 
while simultaneously preserving intricate image textures and details. 
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